clear; close all;

% Data from Jesper:
para=zeros(25,6);
para2323=[0.23582493 2.09349355 0.15563716 1.560913   0.60926898 0.51701385];para(1,:)=para2323;
para4141=[0.18049477 2.92908224 0.49502099 0.99695761 0.3227105  0.36116294];para(7,:)=para4141;
para4242=[0.48796853 0.98194874 0.31579653 0.35741207 0.18798549 2.71874806];para(13,:)=para4242;
para7373=[0.57289615 0.5763908  0.32394661 2.46358374 0.0739523  0.27722376];para(19,:)=para7373;
para7474=[0.32825 2.54931 0.09219 0.29182 0.5811 0.59231];para(25,:)=para7474;
para2341=[0.5026924  1.22895263 0.13951036 7.50424076 0.41849001 0.42091243];para(2,:)=para2341;
para2342=[0.18396472 8.70291804 0.49897344 1.25387083 0.42609115 0.4247833 ];para(3,:)=para2342;
para2373=[0.47408861 0.92079566 0.19860103 3.19599153 0.32249351 0.40575131];para(4,:)=para2373;
para2374=[0.47543791 0.94241702 0.19384435 3.38236532 0.33063109 0.40745887];para(5,:)=para2374;
para4142=[0.19020878 2.65056725 0.32076908 0.36142452 0.47789479 0.98036749];para(8,:)=para4142;
para4173=[0.50494081 0.91672469 0.29923796 0.3728018  0.1826209  3.21613753];para(9,:)=para4173;
para4174=[0.5085639  0.91518537 0.29670954 0.37129392 0.18666422 3.28455076];para(10,:)=para4174;
para4273=[0.51344146 0.98119957 0.16477545 3.69293833 0.31816071 0.37914886];para(14,:)=para4273;
para4274=[0.51684534 0.96288639 0.30904771 0.37485806 0.16779276 3.59037981];para(15,:)=para4274;
para7374=[0.07833469 0.27911064 0.3294006  2.4745664  0.58841741 0.58199911];para(20,:)=para7374;

% output to file to be copied to GPUMD code:
fid=fopen('zbl_para.txt','w');
for n1=0:4
    for n2=0:4
        index=min(n1*5+n2+1,n2*5+n1+1);
        fprintf(fid,'{%.8ff, %.8ff, %.8ff ,%.8ff, %.8ff, %.8ff},\n',para(index,:));
    end
end
fclose(fid);

